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We discuss the set of constraints for Einstein-aether theories, comparing the flat background case 
with what is expected when the gravitational fields are dynamic. We note potential pathologies 
occurring in the weak gravitational field limit for some of the Einstein-aether theories. 
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I. INTRODUCTION 



There has been significant interest in recent years in gravitational field theories that involve a fundamental vector 
field as well as the usual metric field. These include the paramerized family of Einstein- aether theories [l| and the 
TeVeS theory [2j . In addition, the Horava theory [3| can be viewed as a particular case of Einstein-aether theory, so 
it too is effectively a gravitational field theory with a fundamental vector field coupled to a metric field. A number of 
years ago, Isenberg and Nester [J] noted a possible difficulty for certain field theories involving a vector field coupled 
to a spacetime metric field. 1 This difficulty involves the number of constraint equations the theory imposes on the 
choice of initial data for the theory, and how this number changes if one compares (i) the field theory with both the 
1 vector field and the metric field fully dynamic; and (ii) the field theory with the spacetime metric fixed and flat, and 
, only the vector field fully dynamic. 

■ To understand the problem, it is useful to first recall the situation for the Einstein-Maxwell field theory. On a fixed 
flat spacetime background, with no dynamical gravitational fields, the Maxwell theory imposes the two constraints 
D C E C — and D C B C = on the choice of the Maxwell initial data {E a , B a }. For the fully dynamic Einstein-Maxwell 
theory, there are 2 + 4 constraints — including the four Einstein constraints R — K cd K c d + K 2 = rj[E c E c + B C B C ] and 
J> ■ D C K^ — D a K = r/[i? x B]„ as well as the two Maxwell divergence constraints listed above — restricting the choice 
of the Einstein-Maxwell data {h a b, K a b, E a , B a }. 2 This (Einstein-Maxwell theory) is the normal situation for a well- 
behaved theory: there are I = 2 constraints for the fixed background (Maxwell) theory, there are I + 4 constraints in 
the dynamical gravitational field (Einstein-Maxwell) theory, and the constraints behave well in the weak gravitational 
field limit. 

If alternatively one considers the Einstein-vector field theory with the spacetime action principle 3 

S[U a ,g a0 ] = J d A x^g-[K - p[-^W a U^^U a - ^m 2 U a U a ] }, (1) 
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1 The same problem can arise for field theories involving tensor fields coupled to a metric field. 

2 Here E a and B a are the electric and magnetic spatial vector fields, h a j, is the spatial metric field, K a i, is the extrinsic curvature, K 
is the trace of K a i, (with respect to h a i,), D a is the covariant derivative compatible with h a ^, R is the spatial scalar curvature, n is a 
coupling constant, and lower case Latin indices range over the three spatial directions. 

3 Here U a is a spacetime vector field, g a p is the spacetime metric field with the compatible spacetime covariant derivative V a and 
spacetime scalar curvature 1Z, p is a coupling constant, m is a ("mass") constant, and lower case Greek indices range over four spacetime 
directions. 
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one finds the following alternative situation (See [4j for details): For the flat background theory, if one does a space 
+ time decomposition U a — > {U^, U a } of the fields, and if one defines the conjugate momenta {IIx, n a }, then there 
are six constraint equations on the initial data, of the following form 

IT = D a U~ L (2) 

D b H ± = D b D c U c + m 2 U b . (3) 

These constraints effectively reduce the theory to one (vector field) degree of freedom. However, for the fully dynamic 
Einstein-vector version of this theory, instead of having 6+4 constraints on the initial data {h a b, K a j„ U , H±, U a , Il a }, 
one finds that there are only four constraints (analogous to, but much messier algebraically than the usual four Einstein 
constraints). There are then effectively four vector field degrees of freedom, plus two gravitational degrees of freedom. 
However, if the Einstein-vector theory is well behaved then its weak field limit must reduce to vacuum linearized 
gravity for the metric and the flat background theory for the vector field. Consequently, the mismatch between the 
number of constraints shows that the Einstein-vector theory behaves poorly as the gravitational field approaches 
flatness. 

There are two features of the Einstein-vector theory with action S[U a , g a p] (as in (JIJ above) which lead to this 
problem with constraints (and correspondingly with degrees of freedom) in the weak field limit: (i) the flat background 
version of the theory has constraints; and (ii) the theory is "derivative-coupled" in the sense that the non-gravitational 
part of the action involves covariant derivatives. This first feature is crucial, since the difficulty of interest involves the 
"loss" of constraints that occurs as the dynamical gravitational field is turned on. The second feature is also crucial, 
since it is the presence of terms involving K a b in the non-gravitational part of the action principle (following its 3 + 1 
decomposition) that leads to the loss of constraints. We note that the Einstein-Maxwell theory has constraints in 
its fiat background version, but is not derivative-coupled; consequently it avoids this problem. The Einstein-massive 
Klein-Gordon theory, whose action is the same as S[U a , g a p] above except that the term —^\7 pU a \7^U a replaces the 
term — a UpV^U a in (JTJ, is derivative-coupled, but has no constraints in its flat background version; consequently 
it too avoids the problem. 

Where do these constraints in the flat space theory come from? The easiest way to see this is to consider the process 
of finding a Hamiltonian formulation of the theory. Recall that to find a Hamiltonian formulation of a theory, one 
begins with a Lagrangian and varies the Lagrangian with respect to the velocities to obtain the momenta. The next 
step is to invert the velocity-momentum relation to obtain the velocity in terms of the momentum. However, that 
step can fail if the velocity- momentum relation is not invertible. Such a lack of invertibility gives rise to constraints, 
and it is precisely this situation that gives rise to the constraints found in the flat space theories studied in [ij . 

What about the Einstein-aether theories? As we discuss below in Section [III these theories form a four-parameter 
set, with (essentially) every theory in the set involving derivative-coupling. Hence, the first step in applying the 
methods of [i[ to the Einstein-aether theories is to determine which of the theories, in their flat background versions, 
involve constraint equations. We do this in Section IIII1 working with the theories in Hamiltonian form. We find 
that some of the flat background theories involve no constraints ("safe" theories), some of them involve constraints 
regardless of the initial data ( "endangered" theories) and the rest involve constraints for some ranges of data, but not 
for others ("conditionally endangered" theories). In Section ITVl we consider the implications of the results of Section 
IIIII for Einstein-aether theories. This is not as straightforward as for the vector theories of [4| because the weak field 
limit of Einstein-aether theories is not the flat background theory, but rather a theory that still couples the vector field 
to the metric perturbation. Nonetheless we find that the theories that we call "endangered" really are pathological 
in the weak field limit. We also consider the possible implications of the "conditionally endangered" theories. 

II. EINSTEIN-AETHER THEORIES 

The Einstein-aether theories make up a four-parameter family of classical metric-vector theories specified by the 
spacetime action principles 

S Ea [u a ,g a p] = J ' d i x^\n-K^V a u^V fi u v + A(l + «? Q(3U f V 3 )}, (4) 

where our conventions on the spacetime metric and indices are as above, where u a denotes a spacetime vector field, 
where A is a Lagrange multiplier, and where the aether-parameter matrix takes the form 

K$ = ci3 Q/ V + <»6«65 + c 3 d^ - c A u a uPg^. (5) 
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with (constant) parameters {01,02,03,04}. We note that if 02 = 03 = 04 = and A = 0, then this is the Einstein- 
Klein-Gordon (massless) theory; if ci = c 2 = c 4 = and A = 0, then we have essentially the Einstein-vector theory 
([T]) discussed above. If c\ + C3 = C2 = C4 = and A does not vanish, then this is the particular version of the 
Einstein-aether theory discussed in [5|. 

The presence of the Lagrange multiplier term serves to enforce the a priori restriction 

g a pu a uP = -1, (6) 

which says that the vector field must always be a unit-length timelike vector field. While this restriction is crucial to 
the physical application of the Einstein-aether theories, it does not appear to have a major qualitative effect on the 
issue under discussion here. 

As noted above, our focus here is on the Einstein-aether theories with a fixed flat background. In that case, with 
Vap representing the Minkowski metric, the action reduces to 

S E a {f iat)[u a } = J d 4 x{-K$V a u^ p u" + \(l+ VaP u a uf 3 )Y (7) 
the aether-parameter matrix takes the form 

K Uiat)^ = ci?rV + t%SZ6g + c 3 52d? - C4U a u% v , (8) 
and the unit length condition (J6j) can be written as 

ri aP u a u fi = -1. (9) 

One way to determine the presence of constraint equations in the initial value formulation of the (flat background) 
aether theories is to first obtain the spacetime covariant field equations for u a (by varying the action S^aifiat) 
then rewrite these equations in 3 + 1 form, and finally search among the 3 + 1 equations for explicit constraints. The 
Hamilton-Dirac approach provides a more systematic way to find the constraints. Since we are primarily concerned 
here with determining for which of the four-parameter family of aether theories there are constraints (rather than 
determining the exact nature of these constraints) , in fact we need only carry out the first part of the Hamilton-Dirac 
analysis, as we discuss in the next section. For completeness, we summarize the covariant approach in the Appendix. 



III. HAMILTON-DIRAC ANALYSIS 



The first step in the Hamilton-Dirac analysis of a given field theory is to carry out a 3 + 1 (space + time) decom- 
position of the fields and their derivatives, and to substitute the resulting expressions into the Lagrangian for the 
given theory. For the present case, we choose a standard slicing for the flat background spacetime with n a as the unit 
time-like normals, and introduce the tensor 

h a/ 3 := g a p + rianp . (10) 

We note that spacetime tensors which are orthogonal to n a can also be considered as spatial tensors. Viewed in this 
way, h a p (which we can also write as h a b) is both the spatial metric and a projection operator that takes spacetime 
tensors to their spatial part. Using projection with respect to h a p and contraction with n a we can decompose any 
spacetime tensor into spatial tensors. In particular, the four vector u a can be decomposed as 

u a = Vn a +w a , (11) 

where V is a scalar and w a (which we can also write as w a ) is a spatial vector. The unit-length condition @ then 
implies that V is not an independent field but is instead determined by w a through the relation 

V 2 = l + w a w a . (12) 

Similarly, we can decompose the spacetime derivative of u a as 

V a it£ = ~n a npV + npD a V — n a wp + D a wp , (13) 

where an overdot denotes the derivative with respect to time, and D a (which we can also write as D a ) denotes the 
spatial covariant derivative operator. 
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Using eqn. (|T3| in eqn. we find that the Lagrangian corresponding to the action Ssafflat) [ ua ] takes the form 

L = Mw a w a - (M + c 23 )V 2 + 2w a (-c 3 D a V + c 4 Vw b D b w a ) 

-2V(c 2 D a w a + c A Vw a D a V) + Z . (14) 

Here the quantities M and Z are given by 

M := Cl +c 4 U 2 , (15) 
Z := ci(£> a y£>°V - D a w b D a w b ) - c 2 {D a w a f - c 3 D a w b D b w a 

+c 4 ((w a D a w b )(w c D c w b ) - (w a D a V) 2 ), (16) 

and we are using the notation of [l| in which c 23 is an abbreviation for c 2 + c 3 (with corresponding abbreviations for 
any other sums of the Cj). 

While the condition (fT2"]) is a constraint, it is one which (if we restrict attention to non- negative V) can be eliminated 
algebraically. More specifically, we adopt the point of view that (flat background) aether theory is a theory of the 
spatial vector field w a , and that any terms in the Lagrangian that depend on V are to be viewed as simply more 
complicated functions of w a given by substituting y/1 + w a w a for each occurrence of V. In this way, the aether theory 
Lagrangian has no a priori constraints and correspondingly no need for the Lagrange multiplier term. In particular 
we may substitute 

V = V^waw 11 (17) 



along with V = y/1 + w a w a into the Lagrangian expression (|14|l , thereby obtaining 

L = q ab w a w b + 2w a B a + Z, (18) 

where the tensor q ab and vector B a are given by 

q ab := Mh ab - (M + c 23 )V~ 2 w a w b , (19) 
B a := c A V 2 w b D b {V- 1 w a ) - c 3 D a V - c 2 V- 1 w a D b w b . (20) 

This is the 3 + 1 form of the Lagrangian expression for the aether theories which we work with here. 

The next step of the Hamilton-Dirac analysis is to calculate the expressions for the momenta p a conjugate to the 
fields w a as functions of w a and w a , and then attempt to invert these expressions to obtain new expressions for w a 
as functions of w a and p a . It is at this step — a key step in the specification of the Legendre transform which maps 
from a Lagrangian formulation to a Hamiltonian formulation of a given theory — that one finds out if the theory has 
any constraints on the choice of initial data sets. Such constraints exist if and only if one cannot invert the map from 
p a (w,w) to u> a (u>,p). 

To determine the full explicit set of constraints for a given theory via the Hamilton-Dirac analysis, one proceeds 
from this step to construct the Hamiltonian (or set of Hamiltonians) for the theory, incorporating the constraints 
obtained from the non-invertibility just discussed, and one calculates the time derivatives of the constraints using 
this Hamiltonian. We leave the remaining details to other references (see, e.g., @). Here, our only real concern is to 
determine the sets of choice of the Einstein-aether parameters {ci, c 2 , c 3 , C4} for which there are constraints and those 
choices for which there are not. 

Calculating the conjugate momentum p a by varying the Lagrangian (|18[) with respect to w a , 

Pa Sw a =SL = Sw a 2(q ab w b + B a ), (21) 

we obtain 

Pa = 2(q ab W b + B a ). (22) 

From eqn. (|22[) . one immediately sees that the momentum- velocity relation can be inverted if and only if the matrix 
q ab is invertible. It is easy to see that if q ab is not invertible then there is a constraint, since for any vector s a for 
which q ab s b = it follows from eqn. (|2"2l that 

s a (p a - 2Ba) = 0. (23) 

To find the conditions under which q ab is invertible, we note (from eqn. (1191) ). that any vector orthogonal to w a is 
necessarily an eigenvector of q ab , with eigenvalue M. Furthermore, it also follows from eqn. (|19[) that w a is generally 
an eigenvector of q ab with eigenvalue V~ 2 N, where 

N = M - c 23 w a w a . (24) 
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Therefore q ab is invertible if and only if MN =^ 0. This result immediately gives rise to a classification of (flat 
background) aether theories, depending on whether MN vanishes at all points of configuration space ("endangered" 
theories), at no points of configuration space ("safe" theories), or at only some points of configuration space ("condi- 
tionally endangered" theories). We would like to express this classification directly in terms of the parameters of the 
theory. From cqns. (fT5|) and (|24[) we determine that 



MN = (ci4 + c 4 w a w a )(ci4 + (c 4 - c 23 )w a w a ). (25) 
Consequently, we find that the endangered theories satisfy either the condition 

ci = c 4 = 0, (26) 

or the condition 

ci ^ 0, c 4 = c 23 = -ci. (27) 

In particular, the theories studied in Q are analogs for non-unit massive fields of endangered aether theories. The 
safe theories satisfy the conditions 

W0 , *->0, ^^£>0. (28) 

C14 C14 

Any theory that is neither endangered nor safe is conditionally endangered. This completes our classification. 

For the safe theories, it is straightforward to proceed to construct a Hamiltonian. In particular, so long as MN ^ 0, 
we calculate from eqn. (TT91 that 

(TT 6 = J^( Nhab + ( M + W). (29) 
It then follows that the velocity-momentum relation can be inverted to yield 

w a = ±(q- 1 ) ab ( Pb -2B b ). (30) 
Using the non-constrained definition of the Hamiltonian, we have 

H = Pa w a - L = i(<rTV - 2B a )( Pb - 2B b ) - Z. (31) 

Then using eqn. (f29|) . we obtain 

H = 1^( Nhab + (M + c 23 )w a w b )( Pa - 2B a )( Pb - 2B b ) - Z. (32) 

A Hamiltonian which incorporates the collateral constraints can also be constructed for the endangered and con- 
ditionally endangered theories; the details of this construction (peripheral to our concerns here) are discussed in 



IV. CONCLUSIONS 

What are the possible physical implications of our classification of (flat background) aether theories? In Q it is 
shown that the extra constraints in the flat space vector theories studied in that paper are not present in the gravitating 
versions of those theories, as a consequence of the presence of many terms containing h ab appearing throughout the 
3 + 1 expression for the Lagrangian of the coupled theories. We have not performed the corresponding analysis for 
the Einstein-aether theories, endangered or otherwise. Their Lagrangians do, however, contain a similar number of 
h ab terms. We note in addition that it has been shown that the diffeomorphism invariance of Einstein-aether theories 
gives rise in each of them to at least four constraints, just as in general relativity. 0, HI- We are consequently led to 
believe that in general the Einstein-aether theories have four (and no more) constraints. However, unlike the case 
of [H even if we were to show that the Einstein-aether theories have exactly four constraints, that would not allow 
us to conclude that the weak field limit of the endangered theories is pathological, because the weak field limit of 
Einstein-aether theories is not the non-gravitating aether theory, but still contains a residual coupling between the 
vector field and the metric perturbation. Nonetheless, the presence of constraints in the flat space endangered theories 
at least leads to the suspicion that there is something wrong with the corresponding Einstein-aether theories. And 
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that suspicion turns out to be justified. The endangered theories all have c%4 = 0. (see eqns. (j2"6]) - (f27|) ). However, as 
noted in [10( Einstein-aether theories with C14 = have infinite speed for the spin and spin 1 modes of the theory. 
Thus the endangered flat space theories serve as a diagnostic of a pathology in the corresponding gravitating theory. 

What about the conditionally endangered theories? Here the number of constraints changes at particular "bad" 
points of configuration space. Certainly, this should lead one to worry about the physical viability of the flat space 
theory itself, since even if one were to restrict initial data to a "good" region of configuration space (where no 
constraints are present) there remains the possibility of the dynamics causing the system to evolve to a "bad" part 
of configuration space, thus resulting in an ill-defined or singular evolution. However, it is not clear that these 
properties of the flat space theory give any cause to worry about the gravitating theory: note that the bad points of 
configuration space typically have w a ~ 1 and are thus far from the weakly gravitating case. Thus the pathologies 
of the conditionally endangered theories cannot be considered as a reliable guide to the behavior of the Einstein- 
aether theory with the same Cj. Nonetheless, a non-gravitating theory with singular evolution may give rise to naked 
singularities when coupled to gravity. This raises the question of whether gravitational collapse gives rise to naked 
singularities in some of the Einstein-aether theories. This question could be addressed numerically using e.g. the 
methods of [ll| 
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Appendix A: Covariant field equations 

Varying the Lagrangian of eqn. (J5J) with respect to the aether field u a yields the field equation 

V Q J a /3 + Au ( 3 + c 4 a ct V ; 3U Q = 0, (Al) 
where J a p = K apj p v V pU v and ap — u a \7 a up. Contracting eqn. (|A1[) with u' 3 and solving for A yields 

A = /V„ J a p + c 4 a a a a (A2) 

and substituting this into eqn (|Aip yields 

= V a J a p + c 4 a a \7 p u a + u^{u'V a J a 1 + c 4 a a a a ). (A3) 

However, we have 

J a p = Cl V a u + c 2( 5 q /3 V 7 m 7 + c 3 V p u a - c 4 u a a p . (A4) 

So after some straightforward but tedious algebra the field equation becomes 

= Cl [\7 a \7 a u p - m /3 V q u 7 V q m 7 ] + c 23 [V ( 3V 7 u 7 + ii (3 ii 7 V 7 V Q M Q ] 

+C4[a a V pu a + 2u/3d a a a — u a ~S/ a ai3 — apV a u a ]. (A5) 

We now want to produce an evolution equation for w a by substituting the decomposition of u a of eqn. into 
eqn. (|A5j) . At first it might then seem that the field equations might be an overdetermined system, since they seem to 
provide four equations for the three components of w a . However, the field equation vector is automatically orthogonal 
to u a so it suffices to impose the spatial projection of the field equations since spatial part of the field equations along 
with orthogonality to u a implies time part of the field equations. Some straightforward but tedious calculation then 
shows that the spatial projection of eqn (|A5|) becomes 

Mw b -c 23 VVw b = R b , (A6) 
where M is given by eqn. (|15p and R b is given by 

R b = Cl [D a D a w b - w h (V 2 - D a VD a V - w a w a + D a w c D a w c )] 
+c 23 [D b V + D b D a w a + w b (VD a w a + w a D a (V + D c w c ))] 
+c 4 [-(F + D a w a )L b - VVw h - Vw a D a w b - Vw a D a w b 
-w a D a L b - (VV + w a D a V)D b V + L a D b w a + 2w b (L a L a - (VV + w a D a vf )]. (A7) 
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Here the vector L a is the spatial projection of the acceleration vector a a . A straightforward computation yields 

L a = Vw a + w b D b w a . (A8) 

The field equation is not quite in the form that we would like, namely with w a alone on one side of the equation 
and no quantities with second time derivatives on the other side of the equation. However, we can easily put it in 
that form by using eqn. (1121) and its time derivatives. In particular, the first time derivative of eqn. (|12[) yields eqn. 
(|17|) while a second time derivative yields 

V = V- 1 [w a w a + w a w a - V 2 }. (A9) 
Then contracting eqn (|A6[) with w a and using eqn. (|A9|) we find 

NVV = R a w a + M{w a w a - V 2 ), (A10) 
where the quantity N is given by eqn. (|24j) . Then on substituting this result into eqn (|A6[) we obtain 

MNw a = NR a + c 23 [R b w b + M(w b w b - V 2 )]w a . (All) 
We note that if MN ^ 0, eqn. (|A11|) is an equation of motion for w a , while if M N = it is a constraint equation. 
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